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ABSTRACT 

Aims. Polycyclic Aromatic Hydrocarbons (PAHs) have been detected toward molecular clouds and some young stars with 
disks, but have not yet been associated with embedded young stars. We present a sensitive mid-infrared spectroscopic 
survey of PAH features toward a sample of low-mass embedded young stellar objects (YSOs). The aim is to put 
constraints on the PAH abundance in the embedded phase of star formation using radiative transfer modeling. 
Methods. VLT-ISAAC L-band spectra for 39 sources and Spitzer IRS spectra for 53 sources are presented. Line intensities 
are compared to recent surveys of Herbig Ae/Be and T Tauri stars. The radiative transfer codes RADMC and RADICAL 
are used to model the PAH emission from embedded YSOs consisting of a pre-main-sequence star with a circumstellar 
disk embedded in an envelope. The dependence of the PAH feature on PAH abundance, stellar radiation field, inclination 
and the extinction by the surrounding envelope is studied. 

Results. The 3.3 fim PAH feature is undetected for the majority of the sample (97%), with typical upper limits of 5 
xl0~ 16 Wm" 2 . One source originally classified as class I, IRS 48, shows a strong 3.3 [im feature from a disk. Compact 
11.2 nm PAH emission is seen directly towards 1 out of the 53 Spitzer Short-High spectra, for a source that is borderline 
embedded. For all 12 sources with both VLT and Spitzer spectra, no PAH features are detected in either. In total, PAH 
features are detected toward at most 1 out of 63 (candidate) embedded protostars (< 2%), even lower than observed 
for class II T Tauri stars with disks (11-14%). Models predict the 7.7 /j,m feature as the best tracer of PAH emission, 
while the 3.3 fim feature is relatively weak. Assuming typical class I stellar and envelope parameters, the absence of 
PAHs emission is most likely explained by the absence of emitting carriers through a PAH abundance at least an order 
of magnitude lower than in molecular clouds but similar to that found in disks. Thus, most PAHs likely enter the 
protoplanetary disks frozen out in icy layers on dust grains and/or in coagulated form. 

Key words. Stars: pre-main sequence - planetary systems: protoplanetary disks - Circumstellar matter - Astrochemistry 



1. Introduction 
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^P* Polycyclic Aromatic Hydrocarbons (PAHs) have been ob- 
^ served toward a wide range of astrophysical en vironments 

k> ( |Allamandola et aT| 1989 Pceters et al. 20021, including 
the interstellar medium (ISM) and star-forming regions, 
first hinted at by the discovery of widespread broad emis- 
sion features in the mid-infrared. Since PAHs are a good 
tracer of UV radiation, they are an indirect probe of star 
formation in the high opacity environments of molecular 
clouds as well as circumstellar disks. Recent spectroscopic 
studies have detected PAHs toward a significant fraction, 
11-14% and 54% of low-mass and inter mediate mass pre- 
main-sequence (PMS) stars respectively (Acke & van den 



Ancker |2004| |Geers|[2~007[ chapter 2), and recent ground- 
based high spatial resolution observations show evidence for 
t he PAH emission to originate from the circumstellar disks 



(Rcs sler fc Barsony||2003| |Habart et aT.||2006| |Geers et al. 
2007 a|b[ ). PAHs are also prominently seen toward photon 
dominated regions (P DRs) in dense clouds ex posed to mas- 
sive young stars (e.g., |Verstraete et al.| |l996). However, to 
this date, no significant PAH emission has been reported to 
be directly associated with the earlier class 0-1 phase, when 
a (proto)star with circumstellar disk is still embedded in 
an envelope of gas and dust. Deeply embedded high-mass 
young st ellar objects (YSO's) also lack pr ominent PAH 
( van Dishoeck fc van der Tak||2000[ ). 
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Based on observations obtained at the European Southern 
Obse rvatory, Paranal, Chile, within the observing pro gram 



164.1-0605 (VLT-ISAAC). 



PAHs are believed to form in the outflows of carbon- 
rich Asymptotic Giant Branch (AGB) stars, which deposi t 
them in the interstellar medium (e.g., jHelling et al.| l996). 
Their ubiquitous presence in the ISM shows that at least 
larger PAH molecules (of 50-100 carbon atoms or more) 
will survive this phase, until being incorporated in molec- 
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ular clouds. The amou nt of carbon locked up in PAHs is 
" 5 dHabart et al.[ 



[C/H]p AH 5 x 10- 
PAH abundances of 



20041, corresponding 
ative to H, for PAH 



5 x 10"' re 

molecules of 100 carbon atoms, making them among the 
most abundant molecules after H2 and CO. 

PAHs can play an important role in the star-forming en- 
vironments. As large molecules, they provide efficient heat- 
ing of the gas in both the ISM and circumstellar disks. 
Small dust particles like PAHs have also been suggested as 
the main formation site of molecules such as H2 and water 
in the later evolutionary phases of circumst ellar disks when 
classi cal grains have grown to large sizes (Jonkheid et al. 
20061. 



The shape and relative strength of the PAH features in 
the 6-9 /Ltm region has been shown to vary between var- 
ious astrophysical environments from AGBs, the ISM to 
the class II disk sources. This has provided evidence that 
the emission characteristics of PAHs are sensitive to local 
physical conditions and that in terstellar PAHs und ergo pro- 
cessing in space environments ( Peeters et al.|2 002). Further 
support for this comes from recent values for the PAH abun- 
dance in disk surface layers, which is a factor of 10-100 
lower than in the IS M for sources with sufficient UV radi- 
ation to excite them ( Geers et al. |2006[ | . PAHs frozen out 
on icy grains can chemically react with other molecules to 
form a large variety of complex species, including pre-biotic 



molecules ( Bernstein et al.| 1999] Ehrenfreund et al.||2006 l. 
A study of PAHs toward embedded low-mass young stars 
is important to find out what happens to PAHs and how 
large a role they may play in these embedded objects. 

In this paper we present the first mid-infrared spectro- 
scopic survey for PAH emission from embedded low-mass 
young stellar objects using the Infrared Spectrometer And 
Array Camera (ISAAC) mounted on the ESO Very Large 
Telescope (VLT) and the Infrared Spectrograph (IRS) on 
board the NASA Spitzer Telescope, in a wide range of 
nearby star-forming regions. We subsequently use a radia- 
tive transfer code to model PAH emission from young stars 
with disks embedded in an envelope, and discuss the ob- 
served detection rate in the context of PAH abundance, 
the strength of UV emission from the protostar and the 
extinction from the envelope. 



2. Observations and data reduction 

A large sample of low-mass embedded objects was se- 
lected from two previous mid-infrared studies. A set of 
39 sour ces was selected from a VL T-ISAAC 3-5 /an band 
survey 



( van Dishoeck e t al. 2003 1, of which 32 were pre- 
viously presented in Pontoppidan et ah] ( 2003 1 and three 



additional sources in Thi et al. ( |2006p ~T i br this survey L- 
band spectroscopy was obtained with ISAAC, the Infrared 
Spectrometer And Array Camera, installed at the VLT 
Antu (UT1) at ESO's Paranal Observatory in Chile, in the 
low resolution (R = A/AA = 600-1200) spectroscopic mode 
in the spectral domain 2.8-4.2 /zm using a 0.6" x 120" slit. 
The telescope was operated using a standard chop-throw 
scheme with typical chop-throws of 10-20". A full descrip- 
tion of observation and reduction techniques is given in 
Pontoppidan et al. (20031. These sources were required to 
have a rising Spectral Energy Distribution (SED) in the 
mid-infrared as well as the presence of an H2O ice feature 
at 3 /Ltm. The selected sources are listed in Table [AT2| 



In addition, a sample of 53 Spitzer Space Telescope 
Infrared Spectrograph (Spitzer IRS) Short-High (10-20 /Ltm, 
R -600) spectra and 33 Short-Low (5-14.5 /Ltm, R -100) of 
class I sources were obtained in the context of the "Cores 
to Disks" (c2d) Legacy program ( Evans et al.|2003 ). These 
sources were selected for showing the silicate 10 /im fea- 
ture in absorption, and the list includes a small number 
of known or candidate edge-on disks, which are labelled in 
Table A. 2 Most of them have spectral slopes a between 2 



and 24 /im > 0. The vast majority of class and I sources in 
the c2d sample were fo und to have luminosities of Lboi < 3 
L ( Evans et al. 2008[ ) indicative of low mass young stars 
and consistent with the observational fact that the major- 
ity of the young s tars in the regions studied by c2d a re K- 
and M-type stars (|Luhman|20"07l |Qliveira et al.|2008| ). The 
reduction of the Spitzer spectroscopy was performed us- 
ing the optimal PSF extraction technique developed within 
the c2d program, to allow separation of the compact sourc e 
and extended emission components (Lahuis et al. 2007). 
The source size is determined from the width of the PSF 
function fitted to the source, compared to the width of the 
PSF function fit for standard calibrator stars. Comparison 
of the optimal PSF extraction and the full aperture extrac- 
tion provides a direct estimate of any potential extended 
emission. An application of the same technique to remove 
backgroun d PAH features ha s been given for the source 
( |2006[ see their Fig. 4). The details 
reduction procedures are described 



VSSG 1 in Geers et al 



of the observations anc 



Lahuis et al. (2006) and Lahuis (2007, thesis Chap. 3). 
In total, the sample contains 80 sources, for which 12 
sources were observed with both ISAAC and Spitzer. Based 
on new 2MASS-Spitzer-MIPS 2-24 /tni spectral slopes 
and/or sub-mm line and continuum data, as well as C O 4.7 
/im spectroscopy, the classification has been refined 



van 



I Kempen et~al . 2008; Enoch et al.|[2008| |van Kem pcn 2008; 
Jorg ensen et al.|2008[ Pontoppidan priv. comm.). These up- 
dates have confirmed 53 embedded sources and 17 sources 
to be disk sources, while for 10 sources the classification 



remains uncertain (Table A. 2) 



3. Results and discussion 

3.1. VLT-ISAAC spectra 

The VLT-ISAAC L-band spectra are presented in Fig. [I] 
Only 1 source in this sample shows a clear 3.3 /im PAH 
emission feature, IRS 48. This source is known to have a 
spatially resol ved disk surrounded by a very tenuous enve- 
lope at most (Geers et al. 2007a I. No clear PAH feature 



is seen toward any of the other sources. To derive upper 
limits, a gaussian with a peak flux equal to the 3er noise 
and FWHM = 0.07 /im is generated and integrated to re- 
trieve the line flux of the largest PAH 3.3 /im feature that 
would still fall within the noise. The results are listed in 
Table E21 

The hydrogen Pf S (3.296 /Ltm), Pf 7 (3.739 /im) and Br 
a (4.051 /tm) are detected toward most sources, with the 
exception of IRAS 11068-7717 and Cha IRN. These lines 
are believed to b e produced in the accreti on column onto 
the central star ( Calvet fc Gullbrin g 1998 ), and/or in the 



wind launched close to the star ( Kraus et al. 2008 1 . The 



broad absorption features by water ice (2.8 /im) and the 
3.47 /tm ice feature can be clearly seen toward the majority 
of the sources. 
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3.0 3.2 3.4 3.6 3.8 4.0 
Wavelength (/-tin) 

Fig. 1. ISAAC L-band spectra of low-mass embedded 
sources. The narrow HI lines at 3.3, 3.74 and 4.05 /xm are 
seen in most spectra. 



3.2. Spitzer spectra 

Out of the sample of 53 sources, 50 do not show any PAH 
features. A selection of Spitzer spectra with PAH non- 
detections is presented in Fig. [2] Silicate, water and CO2 
absorption are the predominant features in these spectra. 
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Fig. 2. Spitzer IRS spectra of a sample of low-mass em- 
bedded sources without PAH detections. For clarity, the 
spectra are normalized to their flux at 17 /im and scaled by 
factors of 1, 3 and 8 for CrA IRS7 B, IRAS15398-3359 and 
IRAS 12553-7651 respectively. 



These have been analyze d by Boogert et al. (2008) and 
P ontoppidan et al.| p008 1. 

For one Spitzer source, Oph GY 23, the optimal extrac- 
tion method produces a SH spectrum with the observed 
11.2 /im PAH features, shown in Fig. [3] A small size is 
derived for this source, suggesting that the PAH emission 
is directly associated with the star. The classification of 
this source remains uncertain. Based on Spitzer IRS, IRAC 
and MIPS fluxes, Oph GY 23, previously classified as class 
I, shows a declining SED suggestive of a class II source. 
However, GY 23 is surrounded by reflection nebulosity and 
deep hot CO gas absorption features have been observed 
towards this source (Pontoppidan, priv. comm.), hinting at 
the presence of remnant envelope material. 

For two sources, VSSG 1 and IRS 14 in Ophiuchus, PAH 
features are detected from the extended background emis- 
sion, which is assumed not to be directly associated with a 
central object. 

Overall, no PAHs are detected toward any of the 53 
confirmed embedded sources. PAHs are detected towards 
one of the 17 confirmed disk sources, IRS 48, and one out 
of the 10 sources with an uncertain classification, GY 23. 



4. Radiative transfer model 

Here we address the question of where the PAH emission 
can arise in a disk + envelope system through radiative 
transfer modeling and what the non-detections imply quan- 
titatively about the PAH abundance. 



4.1. Physical structure 

We use the two-di mensional axisymmetric rad iative trans- 
fer code RADMC ( |Dullemond fc Dominik|2004[ ) to calculate 
the temperature structure and scattering source function 
for an embedded YSO using a Monte Carlo technique. This 
code requires stellar parameters, a stellar radiation field, 
a density structure and a set of dust opacities. Scattering 
is assumed to be isotropic. A module to treat the emis- 
sion from quantum-heated PAH molecules and Very Small 
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Table 1. Parameters of the template model. 
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Fig. 3. Spitzer IRS spectra of GY 23, a low-mass borderline 
embedded source with a PAH detection. 



Grains, previously described in Geers et al. (2006 1, has been 
included. To generate image s and spectra, ray-tracing i s 
performed using RADICAL ( pullemond fc Turolla|[2000 |. 
All spectra and SEDs presented here are scaled to a dis- 
tance to the observer of 150 pc. 

The density structure adopted here follows [Crapsi et al 



(20081 and is comprised of three components: the disk, the 
envelope and the outflow cone. The adopted density struc- 
ture for the disk has a power-law dependence along the 
radial coordinate and a gaussian dependence in height, and 
can be expressed as 



Aiisk(r, 9) 



■ exp 



rcos^ 



H(r) 



where 9 is the angle from the axis of symmetry. The varia- 
tion of scale-height, i.e., the flaring of the disk, is described 
in the function H(r) — r ■ Hq/Rq ■ (r/Ro) 2 ^ 7 , correspond- 
ing to the self-irradiated passive disk of |Chiang fc G oldrcich 
( 1997 1. For the disk in our template model, we fix the outer 
radius R = 300 AU, inner radius R la = 0.1 AU, H = 90 
AU, and the disk mass to = 5 x 10 -3 M , see also 

Table □ 

The envelope density follows the theoretical structure 
for a r otating and collapsing spheroid as derived by |Ulrich 
(1976], defined by 



(r,V) = po 



1 + 



cos 6 
cos 0q 



cos 8 
2 cos #n 



■ cos 6( 



where 9q is the solution of the parabolic motion of an 
infalling particle which is given by r/R lot ■ (cos#o — 
cos6')/(cos6'osin 2 6'o) = 1, i? ro t is the centrifugal radius of 
the envelope, and po is the density of the equatorial plane 
at the centrifugal radius. The outer radius of the envelope 
is fixed to 10 000 AU and the centrifugal radius is set to 300 
AU. The value of po is varied to cover a range of envelope 
mass M env between 0.1 and 1.5 M Q , the template model has 
M env = 1.0 M Q . This is somewhat more massive than the 
average Class I source, which has typically M onv < 0.5 M . 
The inner radius of the envelope is determined by the out- 
flow cavity, which crosses the equatorial plane at ~20 AU. 
The streamline outflow is included by setting the density of 
the regions where cos 9q is larger than cos 15° to the same 



Parameter 


value 


Radiation field 


4000 K blackbody 




+ excess UV field 


L« 


1 L Q 


M* 


1.0 M 


M disk 


5 x 10~ 2 M Q 


-Rdisk.in 


0.1 AU 


-^disk,out 


300 AU 


M cnv 


1.0 M Q 


-^cnv,out 


10000 AU 




300 AU 


00 


15° 


PAH abun. 


5xl0~ 7 w.r.t. H 



value of the density of the envelope at the outer radius. 
This results in a funnel-shaped cavity which is conical only 
at large scales, where it presents a semi-aperture of 15°. 
The envelope introduces an optical extinction of Ay = 33 
at i = 45°, for the template model with M onv = 1.0 M . 

For the radiation field, a blackbody spectrum is taken 
with T e g = 4000 K and a luminosity of 1 L for the tem- 
plate model. In addition, given the strong dependence of 
PAHs on UV excitation and observations indicating the 
presence of UV exce ss around actively accreting young 



low-mass PMS stars (Basri & Bertout 



1989 



Hartmann & 



Kenyon 1990), the radiation field f or A < 0.8 pjn is sub- 



stituted by a modified Draine field (van Dishoeck & Black 



1982), scaled up to ensure continuous overlap at A > 
pm. The resulting spectrum is normalized to a luminosity 
of 1 L . For comparison, a separate model is run without 
the modified radiation field. Any radiation field produced 
by a higher mass (A-type) star would be in between these 
two extremes. The excess Draine UV field corresponds to 
an intensity relative to the mean interstellar radiation field, 
Go, of 8 x 10 4 at a radius of 100 AU at the surface of the 
disk. 

4.2. Treatment of dust and PAHs 

For the optical properties of the dust grain population, we 
adopt the se t of dust mixtures and opacities used in |Crapsi| 
|et al. ( 2008 1 : a mixture of 71% silicates with a size dis tribu- 
tion similar to that of ( Weingartner & Draine 20011 with 
a-WD-parameter equal to -2 and turnover point for radii 
larger than 0.3 /im covered by a mantle of ices with a water 
ice abundance of 3.0 x 10~ 4 relative to H2, together with 
29% of car bonaceous grains with a size d istribution similar 
to that of Weingartner & Draine (20011, characterized by 
a turnover point at radius of 4.5 /im, an awD-parameter of 
-2 and a water ice abundance of 1.5 x 10~ 5 relative to H 2 . A 
second dust population is added without the ice layer, for 
regions where the temperature is higher than 90 K, when 
ices will have evaporated. These are included by running 
the Monte Carlo code once with ices in the whole grid to 
calculate the temperature structure, and then using this in- 
formation to replace the dust opacities where the temper- 
ature is higher than 90 K. The ice features include H2O, 
CO2 and CO, but not more complex species which could 
be respo nsible for (part of) th e observed 6.0 and 6.8 /im ice 
features ( |Boogert et al.|[2008 l. 

PAHs are added as a dust species to the disk and the en- 
velope. The PAH emission is calculated for an equal mix of 
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Fig. 4. Model spectra at i = 45° for template parameters 
listed in Table [l] with PAHs (black) and without PAHs 
(red). Major PAH and absorption features are indicated, 
as well as the wavelength coverage of the presented ISAAC 
and Spitzer IRS observations. The inset shows a blow-up of 
the 3.3 /im feature on the red wing of the water absorption 
band. 



neutral and singly i onized C100H24 molecules, adopting the 
Draine & Li (20011 PAH emission model, using the "ther- 
mal continuous" approximation. Multi-photon events are 
included for the PAH excitation, fol lowing the method out- 
lined by |Siebe nmorgc n~et al.| fl!992| ) . We include the opac- 
ities from |Mattioda e t al.| (1 20050 for near-in f rared wave- 
lengths. Model calculations by | Visser et al.| (2007) show 
that PAHs with N c — 100 have longer lifetimes than the 
lifetime of the disk, and we can therefore safely discard the 
possibility of PAH destruction and keep the PAHs at con- 
stant abundance throughout the entire disk. 

In our templ ate model, th e total PAH abundance is the 
same as used in Geers et al. (2006), as a mass fraction of 
the dust of 0.061, which corresponds to a carbon abundance 
of 5xl0 -5 with respect to hydrogen and af c = 100 PAH 
abundance of 5xl0 -7 with respect to hydrogen. This PAH 
abundance is divided into 50% ionized and 50% neutral 



PAHs, the same as in our previous disk-only models ( Geers 
et al.||2006l|2007bp. 



A spectrum of the template model, with and without 
PAHs, is shown in Figure [4] The main features in the spec- 
trum are the PAH features at 3.3, 6.2, 7.7, 8.6, 11.2 and 
12.7 /im, as well as water, CO, silicate and CO2 ice absorp- 
tion features at 3.08, 4.67, 9.7 and 15.2 /im respectively. 
The wavelength coverage of the ISAAC L-band and Spitzer 
IRS spectroscopy is indicated. 

4.3, Modeling results 

The apparent absence of PAH features toward the majority 
of low-mass embedded class I sources could have a number 
of explanations. 

First, PAH molecules are primarily excited by UV pho- 
tons. The precise shape and strength of the radiation field 
inside an embedded object is not well-known but in the 
class I sources the central source has already formed and is 
the main energy source inside the envelope. The presence 
or absence of excess UV will influence the PAH emission 
features. In particular, the high opacity of the envelope at 



UV and optical wavelengths, which provide the main exci- 
tation of PAHs, will constrain the region where PAHs can 
be excited to a small radius. 

Second, even if the infrared PAH features are produced 
in the disk or inner envelope close to the star, the dust in 
the surrounding envelope may provide too high extinction 
in the optical and mid-infrared for the (PAH) emission to 
escape, especially if the envelope mass is (still) relatively 
high compared to the disk. If an outflow cavity is present, 
the inclination at which the object is observed will affect 
the extinction by the envelope. 

Third, the abundance of small PAH molecules in the 
gas phase may be significantly lower due to freeze-out due 
to the low temperatures and high densities in the interior 
of the molecular core. 

In practice, a combination of the above three causes will 
apply simultaneously. To estimate their effects, models are 
run, varying several parameters including radiation field, 
PAH abundance and the mass of the envelope. 

4.3.1. Luminosity of the central source and presence of UV 
excess 

Models with total stellar luminosity varying between 1, 
3 and 6 L are shown in Fig. [5] In addition, a model with- 
out the Draine field for UV excess is shown for — 6 L Q . 
Increasing the luminosity by a factor of 6 increases the line 
flux of the PAH features by the same amount while the 
PAH feature / continuum ratio is reduced by at most 20%. 
Excluding the excess UV field, while preserving the total 
stellar luminosity, decreases the PAH feature to continuum 
ratios by ~ 3. Interestingly, the 3.3 /an feature is affected 
most because it requires higher energies to be excited than 
the other features. For typical luminosities observed toward 
low-mass embedded protostars, PAH features should be de- 
tectable if they have ISM abundances, even if no UV excess 
is present. 

4.3.2. Mass of the envelope 

Models with the mass of the envelope M cnv varying from 
0.1, 0.5, 1.0 to 1.5 M Q are shown in Fig.[6j at i = 45°. In all 
models the PAH features are clearly present. Increasing the 
envelope mass by a factor of 15 results in a decrease of the 
emission of the central source while the sub-mm emission 
and the strength of the absorption features increase. The 
PAH features decrease in peak flux by a factor of 3-4, but 
are in no case extinguished by the continuous extinction 
or silicate abs orption features. Thi s is co nsistent with the 
conclusions of Manske & Henning ( 1999 1 for higher mass 
YSO's. 

In the template models, PAHs are located in both the 
envelope and the disk. To test the influence of the enve- 
lope further, model setups including PAHs only in the disk 
or only in the envelope were performed. A comparison is 
shown in Fig. [7] for M cnv = 1.0 and 0.1 M Q . Including 
PAHs only in the disk results in a decrease of peak flux of 
the PAH features by a factor ~2. In the SED, the absence of 
PAHs in the envelope leads to less absorption of UV emis- 
sion, which appears stronger. The mass of the envelope has 
no significant effect on the strength of the PAH features, al- 
though the underlying continuum changes. Including PAHs 
only in the envelope also results in a typical decrease of the 
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Fig. 5. Model SEDs (top) and blow-up of spectra (bottom) 
at i = 45° for template model parameters and varying 
between 6, 3, and 1 L (black, red, blue respectively), all 
including UV excess. A model with = 6 L Q without 
UV excess is shown in purple. Major PAH and absorption 
features are indicated. The inset shows a blow-up of the 3.3 
/im feature on the red wing of the water absorption band. 



PAH feature peak flux, by about a factor of 1.5. Even if 
the PAHs are located only in the envelope, PAHs should be 
detectable if they are present at ISM abundance. 

4.3.3. Outflow cavity and inclination 

Embedded class I objects are known to have ou tflows (e.g 



Hogerheijde et al. 1998 Arce & Sargent 20061 so that an 



outflow cavity is included in our model envelope. Because of 
the presence of the outflow cavity, the inclination at which 
the object is observed is important, because at near-pole-on 
orientation one observes directly the central source and the 
disk. At larger inclinations, the envelope will be obscuring 
the disk. 

A template model with the template PAH abundance 
and an envelope mass of 1.0 Mq, seen at varying inclina- 
tion angles between 5 and 85° is shown in Fig. [8] At an 
inclination of i = 5° the stellar radiation field (blackbody 
+ scaled Draine field) is directly visible and no absorption 
features arc present. Between inclination of 5° and 25°, i.e., 
down the cavity and through the envelope, the appearance 
of the SEP and the PAH features chan ges rapidly (e.g., 
Whitney et~al]|2003| |Crapsi et al]|2008|. The emission of 



Fig. 6. Model SEDs (top) and blow-up of spectra (bottom) 
with M env = 1.5 (black), 1.0 (red), 0.5 (blue) and 0.1 (pur- 
ple) Mq at i — 45°. The curve for 1.0 Mq was omitted 
from the SED plot for clarity. The inset shows a blow-up 
of the 3.3 jjm PAH feature. 



the central star becomes obscured and the strength of the 
PAH features decreases by about a factor of 2. At increas- 
ing inclination the PAH features become weaker, but still 
dominate the spectrum. Ice absorption features can be seen 
at 3, 4.2, 6 and 15 /im. At 85°, the disk is observed almost 
edge-on and the PAH emitting regions are largely masked 
by the disk itself. The remaining features in the spectrum 
extracted for an infinitely large aperture are due to scat- 
tered emission, originating from higher up in the disk at- 
mosphere. Observed PAH feature strengths will depend on 
the pointing and orientation of the limited aperture slit on 
the embedded source. 

4.3.4. PAH abundance 

Models with the total PAH abundance (both ionized and 
neutral species) varying between the template abundance, 
and factors of 10, 20 and 100 lower, are shown in Fig. 9j 
The envelope mass in these models is 1.0 Mq. Reducing the 
PAH abundance by a factor of 10 has the straightforward 
result of reducing the PAH peak flux by a factor of 3. The 
3.3 /im feature is notably weak in all models, and already 
at a factor of 10 lower abundance not anymore detectable. 
Decreasing the PAH abundance also leads to an increase in 
the UV emission and continuum emission at far-IR wave- 
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Fig. 7. Model SEDs (top) and blow-up of spectra (bottom) 
with PAHs in both envelope and disk (black), only in the 
disk (red) and only in the envelope (blue), at i = 45°, for 
M env = 1.0 (left) and 0.1 M G (right). The inset shows a 
blow-up of the 3.3 /im PAH feature. 



lengths. The strongest feature, 7.7 /im, is visible down to 
an abundance of 1 x 10~ 8 relative to H, i.e., a value 50x 
lower than the ISM. 

4.4. Summary and caveats 

We find in all models that the 7.7 /urn feature is noticeably 
the strongest PAH feature, with the highest feature to con- 
tinuum ratio. This is in sharp contrast with observational 
and modeling results for non-embedded class II stars with 
disks, where the silicate emission feature largely obscures 
the contribution of the 7.7 and 8.6 /im features. There is not 
any readily visible 7.7 /im PAH feature in our observations. 
Increasing the envelope mass and decreasing PAH abun- 
dance, luminosity or UV excess all lead to weaker PAH fea- 
tures. The inclination has the largest influence on the PAH 
feature strength in the range of inclinations where the line- 
of-sight passes from through the envelope to through the 
outflow cavity. Strong PAH features are predicted to oc- 
cur when observing down the outflow cavity. This situation 
can be recognized by a strong UV and optical part of the 
SED, as well as a lack of absorption features. For typical 
values for the luminosity, UV excess, inclination and enve- 
lope mass, derived from observational studies of embedded 
protostars, and assuming an ISM abundance, the PAH fea- 
tures should be detectable. The low detection rate of PAHs 
toward these sources thus suggests that the abundance of 
PAHs is at least 20-50x lower than in the ISM. Decreasing 
luminosity, UV excess or envelope mass compared with the 
template model changes this conclusion to a typical factor 
of 10-20, compared with the ISM abundance. 

One caveat remains with respect to the effect of scatter- 
ing. The radiative transfer code used here assumes all scat- 
tering to be isotropic. This does not allow us to treat the 
effect of preferentially forward scattering of UV radiation 
through the outflow cavity. In the current implementation, 
our models could overestimate the amount of radiation re- 
ceived by the disk as well as overestimate the strength of 
the PAH features. Thus, our inferred abundance limits are 
on the low side. As a second caveat, an external interstellar 
radiation field is not included in these models. Interstellar 
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20 



Fig. 8. Model SEDs (top) and blow-up of spectra (bottom) 
for i = 5, 25, 45, 65 and 85° (black, red, blue, purple, gray, 
respectively). Curves for i = 25 and 65° were omitted from 
the SED plot for clarity. The inset shows a blow-up of the 
3.3 /tm PAH feature. 



UV radiation from outside the embedded object may be- 
come a significant source of excitation for PAHs in the en- 
velope at large radii and would increase the predicted PAH 
feature strength and thus constrain the PAH abundance 
upper limits to smaller values, going in opposite directions. 

4.5. PAH evolution from clouds to disks 

In both class I and II low-mass objects that we have studied, 
we find very low detection rates for PAHs. Comparison with 
model predictions suggest that for typical conditions in 
class I protostars, the absence of exciting radiation and at- 
tenuation by foreground material can be excluded as reason 
for the absence of the PAH features. In addition, the pre- 
dicted feature / continuum ratio is well within the achieved 
sensitivity of the observations. This thus implies an absence 
of carrier as the reason the for absence of features, and the 
PAH abundance in the embedded protostar phase is esti- 
mated here to be at least a factor of 10-20 times lower than 
for typical (number of carbon atoms Nc — 100) PAHs in 
the ISM. This lower abundance is in the same range a s 



those inferred toward T Tauri stars in Geers et al. ( 2006 ) . 



Thus, even though the detection rate in the Class 1 ob- 
jects is even lower than that found toward disks, the in- 
ferred abundances are comparable. Note that the few Class 
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Fig. 9. Model SEDs (top) and blow-up of spectra (bottom) 
with template PAH abundance (black), and factors of 10 
(red), 20 (purple) and 100 lower (blue) at i = 45°. The 
inset shows a blow-up of the 3.3 /xm PAH feature. 



II disks in which we detect PAHs are mostly G-stars. The 
disk sample in Geers et al. (2006) is more biased toward 
higher mass stars. 

What would cause this lower abundance in the embed- 
ded phase, and is this mechanism the same in the class II 
phase? In a cold dense environment, two possibilities for 
lowering the abundance of a species are recognized: coagu- 
lation or dust growth and freeze-out onto larger grains. 

The freeze-out of PAHs onto grains could explain the 
lower PAH abundance in these objects. PAH features in 
absorption have been observed in the ISM, and laboratory 
spectra of PAHs condensed into H 2 rich ices suggest that 
PAH absorption features should also be observable in astro- 
physical ices, such as in dense clouds ( |Bernstein et al.|2 007 
and references therein) . During this phase, if exposed to U V 
radiation, PAHs in water ice can become significantly ion- 
ized, thus promoting ion-mediated rea ctions in these ices 
with implications for astrochemistr y (lEhrenfreund et al. 
2006 Gudipati fc Allamandola|2006[ ). In the class II phase, 
the direct irradiation of the disk by the central star can 
increase the temperature in the surface layers of the disk, 
evaporating the ice and depositing the enclosed PAHs back 
into the gas, thereby increasing the abundance of the PAHs 
again. Recent observations of edge-on disks show evidence 
for the conti nued presence of a reservoir of ice in the inte- 
rior of disks (Pontoppidan et al. 2005 Terada et al. 2007) 



which can be brought to the surface by vertical mixing. If 
the PAHs in the ices react to form larger carbonaceous ma- 
terials, then only a fraction of PAHs might come back off 
the grain when the ice layer evaporates. This could be a pos- 
sible explanation for the low PAH abundance inferred for 
the majority of T Tauri disk surface layers. In a study of the 
i ce features of the sa me c2d sample of embedded sources, 
( |Boogert et al.|[2~008 l show spectra over the 5-10 //mrange 
where absorptions due to PAHs in ices occur. No obvious 
features are detected, but the laboratory spectroscopy is 
not well known and some PAH bands could be hidden in 
the 6 imq. complex, espe cially on the long wa velength side 
of the C2 component in |Boogert et al. (12008} . Assuming a 



typical (but uncertain) oscillator strength for PAHs, up to 
10% of the carbon could be locked up in PAHs in the ices 
(Boogert, priv. comm.). These limits are consistent with 
the amount of freeze-out found in this work. 

PAHs growing and/or agglomerating to typically lOOx 
larger particles could be an alternative explanation of low- 
ering the PAH abundance. This process could take place in 
the high density regions of the disk and envelope, in both 
the class I and II phase, and does not exclude freeze-out 
as an option. Observationally, the presence of broad emis- 
sion plateaus between 6-9, 11-13 and 15-20 /zm has been 
attribute d to large 200-2000 carbon atom PAHs and PAH 



clusters (Allamandola et al. 1989 Van Kerckhoven et al 



2000 ) , which will lose their characteristic PAH emission fea- 



tures, thus additionally lowering the feature / continuum 
of these features. The detection of individual features to- 
ward T Tauri disks shows that this process is not necessar- 
ily dominant, and that either small PAHs are still released 
from the ice in this phase and/or that small PAHs are still 
being created from, e.g., the destruction of larger grains or 
agglomerates through collisions. 

5. Conclusions 

No PAHs are detected toward any of the 53 confirmed em- 
bedded sources in neither our Spitzer nor ISAAC surveys. 
PAHs are detected at 3.3 /jm towards 1 of the 17 confirmed 
disk sources, IRS 48, and at 11.2 /mi for 1 of the 10 sources 
with an uncertain or borderline embedded classification, 
GY 23. For all 12 sources with ISAAC and Spitzer spectra, 
no PAH features are detected in either. In total, PAH fea- 
tures are detected toward at most 1 out of 63 (candidate) 
embedded protostars (1.6%), even lower than observed for 
class II T Tauri stars with disks (11-14%). 

The low PAH detection rate is compared with radiative 
transfer model calculations. It is shown that the effects of 
envelope mass on the absorption of PAH features is not 
enough to obscure PAH features. The 7.7 /im feature is 
predicted by models as the best tracer of PAH emission, 
while the 3.3 /im feature is relatively weak. 

Assuming typical class I stellar and envelope parame- 
ters, the absence of PAH emission is most likely explained 
by the absence of emitting carriers through a much lower 
PAH abundance in the gas, e.g., due to freeze-out of PAHs 
on icy layers on dust grains or agglomeration. Our inferred 
low abundances are similar to those found for disks around 
T Tauri stars. Thus PAHs are expected to be removed from 
the gas already at earlier stages of star- and planet for- 
mation and enter the disks frozen out on grains. Further 
searches for PAH absorption features are needed to con- 
strain the presence of PAHs in icy grains. 
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Table A. 2. Summary of observations. Line flux in W m~ 2 . 



Target 


RA 




Dec 






class 


Line flux 


Comments 




[J2000] 


[J2000] 








3.3 /Jm 






Perseus 


LDN 1448 1RS1 


3 h 25 m 


09 s . 4 


+30° 46' 


21' 


.7 


embedded 




Spitzcr 




LUrN 1448 IN A 


3 zo 


36 s .5 


-f-OU 40 


2 1 


.2 


embedded 




Spitzcr 




IKAb 0324o+oUU2 


3 h 27 m 


39 s .0 


+30° 12 


59 


.4 


embedded 




Spitzcr 




L1455 SMM1 


3 h 27 m 


43 s . 2 


+30° 12' 


28' 


.8 


embedded 




Spitzcr 




L1455 IRS3 


3 h 28 m 


00 s . 4 


+30° 08' 


01' 


.3 


embedded 




Spitzcr 




IRAS 03254+3050 


3 h 28 m 


34 s . 5 


+31° 00' 


51' 


.1 


embedded 




Spitzcr 




IRAS 03271+3013 


3 h 30 m 


15 s . 2 


+30° 23' 


48' 


.8 


embedded 




Spitzcr 




IRAS 03301+3111 


3 h 33 m 


12 s . 8 


+31° 21' 


24' 


.1 


embedded 




Spitzcr 




Bl-a 


3 h 33 m 


16 s . 7 


+31° 07' 


55' 


.2 


embedded 




Spitzcr 




Bl-c 


3 h 33 m 


17 s . 9 


+31° 09' 


31' 


.0 


embedded 




Spitzcr 




SSTc2d J033327.3+310710 


3 h 33 m 


27 s . 3 


+31° 07' 


10' 


.2 


embedded 




Spitzcr 




HH 211-mm 


3 h 43 m 


56 s . 8 


+32° 00' 


50' 


.4 


embedded 




Spitzcr 




IRAS 03439+3233 


3 h 47 m 


05 s . 4 


+32° 43' 


08' 


.4 


embedded/borderline 




Spitzcr 




IRAS 03445+3242 


gh 47 m 


41 s . 6 


+32° 51' 


43' 


.9 


embedded 




Spitzcr 




Taurus 


LDN 1489 IRS 


04 h 04 m 42 s . 9 


+26° 18' 


56' 


.3 


embedded 


<8.8E-16 


ISAAC 




Orion 


Rcipurth 50 


05 h 40 m 27 s .7 


-07° 27 


32' 


.1 


embedded 


<1.0E-15 


ISAAC 




TPSC 78 


05 h 35 m 14 s . 1 


-05° 23 


38' 


.4 


uncertain 


<2.3E-16 


ISAAC 




TPSC 1 


05 h 35 m 14 s . 5 


-05° 23' 


54' 


.7 


uncertain 


<1.0E-16 


ISAAC 




Vela 


HH 46 / IRAS 08242-5050 


08 h 25 m 43 s .8 


-51° 00' 


35' 


.6 


embedded 


<1.2E-16 


ISAAC, Spitzer 


IRAS 08261-5100 


08 h 27" 


"38 s . 9 


-51° 10' 


37' 


.2 


embedded 




Spitzcr 




LLN 20 


08 h 47' 


n 39 s .4 


-43° 06' 


08' 




embedded 


<1.1E-16 


ISAAC 




LLN 33 


08 h 57" 


"36 s . 8 


-43° 14' 


35' 




embedded 


<1.3E-16 


ISAAC 




LLN 47 


09 h 09 m 25 s . 6 


-45° 22 


51' 




embedded 


<5.3E-16 


ISAAC 




Chamaeleon 


Cod 110 IRS4 


ll h 06 m 46 s .6 


-77° 22 


32' 


.5 


embedded 




Spitzcr 




IRAS 11068-7717 


ll h 08 m 15 s . 1 


-77° 33' 


53' 


.2 


uncertain 


<5.5E-15 


ISAAC 




Cha IRN 


ll h 08 m 38 s .2 


-77° 43' 


51' 


.7 


embedded 


<4.2E-16 


ISAAC 




Cha INa 2 


ll h 09 m 36 s .6 


-76° 33 


39' 




disk 


<5.1E-16 


ISAAC 




IRAS 12553-7651 


12 h 59' 


n 06 s . 6 


-77° 07' 


40' 


.1 


embedded 




Spitzcr 




Ophiuchus 


VSSG1 


16 h 26' 


n 18 s . 8 


-24° 28 


49' 


.5 


disk 


<9.5E-16 


ISAAC, 


Spitzcr 


GSS30 IRS 1 


16 h 26' 


„ 21 b 4 


-24° 23' 


04' 


.2 


embedded 


<8.6E-16 


ISAAC, 


Spitzer 


GY 23 


16 h 26' 


„ ^ 


-24° 24' 


48' 


.2 


embedded/borderline 




Spitzcr 




VLA 1623-2418 


16 h 26' 


n 26 s . 4 


-24° 24' 


30' 


.2 


embedded 




Spitzcr 




IRS 14 


16 h 26' 


n 31 s . 


-24° 31 


05' 


.2 


disk 




Spitzcr 




WL 12 


16 h 26' 


» 44 s . 2 


-24° 34' 


48' 


.4 


embedded 


<2.6E-16 


ISAAC, 


Spitzcr 


OphE-MM3 


16 h 27' 


n 05 s . 9 


-24° 37 


08' 


.0 


disk 




Spitzcr 




GY 224 


16 h 27' 


n 11 s . 2 


-24° 40' 


46' 


.6 


disk 




Spitzcr 




WL 19 


16 h 27' 


n 11 s . 7 


-24° 38 


32' 


.3 


uncertain 




Spitzcr 




WL 20S 


16 h 27' 


n 15 s . 6 


-24° 38' 


45' 


.6 


disk 


<2.0E-16 


ISAAC, 


Spitzer 


WL 20E 


16 h 27' 


n 15 s . 7 


-24° 38' 


59' 


.8 


disk 


<2.1E-16 


ISAAC 




IRS 37 


16 h 27' 


n 17 s . 6 


-24° 28' 


56' 


.6 


embedded 




Spitzcr 




IRS 42 


16 h 27' 


n 21 s . 5 


-24° 41 


43' 


.1 


disk 


<2.0E-15 


ISAAC 




WL 6 


16 h 27' 


n 21 s . 8 


-24° 29 


53' 


.2 


embedded 




Spitzcr 




CRBR 2422.8-3423 


16 h 27' 


n 24 s . 6 


-24° 41 


03' 


.1 


disk 


<8.7E-16 


ISAAC, 


Spitzer 


IRS 43 


16 h 27' 


n 26 s . 9 


-24° 40' 


50' 


.8 


embedded 


<4.4E-15 


ISAAC 




IRS 44 


16 h 27' 


° 28 s . 


-24° 39' 


33' 


.5 


embedded 


<1.3E-15 


ISAAC 




Elias 32 


16 h 27' 


n 2gS 4 


-24° 27 


21' 


.2 


embedded 


<8.4E-16 


ISAAC, 


Spitzcr 


IRS 46 


16 h 27' 


„ 2gs 4 


-24° 39' 


16' 


.2 


disk 


<1.4E-15 


ISAAC, 


Spitzcr 


VSSG 17 


16 h 27' 


n 30 s . 2 


-24° 27 


44' 


.3 


embedded 


<3.9E-16 


ISAAC, 


Spitzer 


IRS 48 


16 h 27' 


„ 37 s 2 


-24° 30' 


35' 


.0 


disk 


3.0E-15 


ISAAC 




IRS 51 


16 h 27' 


n 39 s . 8 


-24° 43 


15' 


.1 


disk 


<9.2E-16 


ISAAC 




IRS 54 


16 h 27' 


n 51 s . 8 


-24° 31 


45' 


.5 


disk 


<3.8E-16 


ISAAC 




IRS 63 


16 h 31' 


n 35 s . 7 


-24° 01' 


29' 


.6 


embedded 


<3.2E-15 


ISAAC, 


Spitzcr 


L1689-IRS5 


16 h 31' 


n 52 s . 1 


-24° 56 


15' 


.4 


uncertain 




Spitzcr 




IRAS 16293-2422B 


16 h 32' 


n 22 s . 6 


-24° 28' 


32' 


.2 


embedded 




Spitzcr 




IRAS 16293-2422 


16 h 32' 


n 22 s . 9 


-24° 28 


36' 


.1 


embedded 




Spitzcr 




RNO 91 


16 h 34" 


n 29 s . 3 


-15° 47' 


01' 


.3 


embedded 




Spitzcr 




Serpens 


SSTc2d J182901. 8+02954 


18 h 29" 


"01 s . 8 


+00° 29' 


54' 


.2 


uncertain 




Spitzcr 




SSTc2d J182916. 2+01822 


18 h 29' 


n 16 s . 2 


+00° 18' 


22' 


.7 


embedded 




Spitzcr 




Serp-S68N 


18 h 29' 


n 48 s . 1 


+01° 16' 


42' 


.6 


embedded 




Spitzcr 




EC69 


18 h 29' 


n 54 s . 4 


+01° 15' 


01' 


.8 


uncertain 




Spitzcr 




SVS 4-2 


18 h 29' 


n 56 s . 6 


+01° 12' 


59' 


.4 


disk 


<6.4E-17 


ISAAC 




Serp-SMM4 


18 h 29' 


n 56 s . 6 


+01° 13' 


15' 


.2 


embedded 




Spitzcr 




EC 82 


18 h 29' 


n 56 s . 9 


+01° 14' 


46' 


.4 


disk 


<4.9E-16 


ISAAC 




EC 90A 


18 h 29" 


n 57 s . 3 


+01° 14' 


03' 


.7 


embedded 


<2.1E-15 


ISAAC 




EC 90B 


18 h 29" 


n 57 s . 3 


+01° 14' 


03' 


.7 


embedded 


<7.4E-16 


ISAAC 




EC88 


18 h 29 m 57 s .6 


+01° 13' 


00' 


.5 


embedded 




Spitzcr 




SVS 4-3 


18 h 29 m 56 s .7 


+01° 12' 


39' 


.0 


disk 


<1.4E-16 


ISAAC 




SVS 4-5 


18 h 29 m 57 s . 6 


+01° 13' 


00' 


.2 


disk 


<1.6E-16 


ISAAC 





Geers et al.: Lack of PAH emission toward low-mass embedded YSOs 



Target 


RA 

[J2000] 




Dec 
[J2000] 




class 


Line flux 
3.3 /im 


Comments 


SVS 4-9 


lg h 2g m 


58 s . 1 


+01° 12' 


39". 4 


disk 


<3.2E-16 


ISAAC 


Scrp-SMM3 


18 h 29 m 


59 s . 2 


+01° 14' 


00". 2 


embedded 




Spitzcr 


Coronis Australis 


R CrA IRS 5A 


19 h 01 m 


48 s . 1 


-36° 57' 


21". 9 


embedded 


<4.2E-16 


ISAAC 


R CrA IRS 5B 


19 h 01 m 


48 s . 1 


-36° 57' 


21". 9 


embedded 


<1.7E-16 


ISAAC 


R CrA IRS 5A+B 


19 h 01 m 


48 s . 


-36° 57' 


21". 6 


embedded 




Spitzcr 


HH 100 IRS 


19 h 01 m 


50 s . 7 


-36° 58' 


09". 6 


embedded 


<8.8E-15 


ISAAC 


R CrA IRS 7A 


19 h 01 m 


55 s .3 


-36° 57' 


22". 


embedded 


<1.5E-16 


ISAAC, Spitzcr 


R CrA IRS 7B 


19 h 01 m 


56 s . 4 


-36° 57' 


28". 1 


embedded/borderline 


<7.1E-17 


ISAAC, Spitzcr 


R CrA IRAS32 


19 h 02 m 


58 s . 7 


-37° 07' 


34". 7 


embedded 




Spitzcr 


Additional sources 


IRAS 13546-3941 


13 h 57 m 


38 s . 9 


-39° 56' 


00". 2 


disk 




Spitzcr 


IRAS 15398-3359 


15 h 43 m 


02 s . 3 


-34° 09' 


06". 8 


embedded 




Spitzcr 


IRAS 23238+7401 


23 h 25 m 


46 s . 7 


+74° 17' 37". 3 


embedded 




Spitzcr 



Geers et al.: Lack of PAH emission toward low-mass embedded YSOs 



13 



References 

Acke, B. & van den Ancker, M. E. 2004, A&A, 426, 151 
Allamandola, L. J., Tielens, G. G. M., & Barker, J. R. 1989, ApJS, 
71, 733 

Arce, H. G. & Sargent, A. I. 2006, ApJ, 646, 1070 
Basri, G. & Bertout, C. 1989, ApJ, 341, 340 

Bernstein, M. P., Sandford, S. A., Allamandola, L. J., ct al. 1999, 

Science, 283, 1135 
Bernstein, M. P., Sandford, S. A., Mattioda, A. L., & Allamandola, 

L. J. 2007, ApJ, 664, 1264 
Boogert, A. C A., Pontoppidan, K. M., Knez, C, et al. 2008, ApJ, 

678, 985 

Calvet, N. & Gullbring, E. 1998, ApJ, 509, 802 

Chiang, E. I. & Goldreich, P. 1997, ApJ, 490, 368 

Crapsi, A., van Dishoeck, E. F., Hogerheijde, M. R., Pontoppidan, 

K. M., & Dullemond, C. P. 2008, A&A, 486, 245 
Draine, B. T. & Li, A. 2001, ApJ, 551, 807 
Dullemond, C. P. & Dominik, C. 2004, A&A, 417, 159 
Dullemond, C. P. & Turolla, R. 2000, A&A, 360, 1187 
Ehrenfreund, P., Rasmusscn, S., Cleaves, J., & Chen, L. 2006, 

Astrobiology, 6, 490 
Enoch, M. L., Evans, II, N. J., Sargent, A. I., et al. 2008, ApJ accepted, 

ArXiv 0805.1075 

Evans, II, N. J., Allen, L. E., Blake, G. A., et al. 2003, PASP, 115, 
965 

Evans, II, N. J., Dunham, M. M., j0rgensen, J. K., ct al. 2008, ApJ 

in press, arXiv:081 1.1059 
Geers, V. C. 2007, PhD in Astronomy, Leiden Observatory - Leiden 

University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
Geers, V. C, Augereau, J.-C, Pontoppidan, K. M., et al. 2006, A&A, 

459, 545 

Geers, V. C, Pontoppidan, K. M., van Dishoeck, E. F., et al. 2007a, 
A&A, 469, L35 

Geers, V. C, van Dishoeck, E. F., Visser, R., et al. 2007b, A&A, 476, 
279 

Gudipati, M. S. & Allamandola, L. J. 2006, ApJ, 638, 286 
Habart, E., Natta, A., & Kriigcl, E. 2004, A&A, 427, 179 
Habart, E., Natta, A., Testi, L., & Carbillet, M. 2006, A&A, 449, 1067 
Hartmann, L. W. & Kenyon, S. J. 1990, ApJ, 349, 190 
Helling, C, Jorgensen, U. G., Plez, B., & Johnson, H. R. 1996, A&A, 
315, 194 

Hogerheijde, M. R., van Dishoeck, E. F., Blake, G. A., & van 

Langevelde, H. J. 1998, ApJ, 502, 315 
Jonkheid, B., Kamp, I., Augereau, J.-C, & van Dishoeck, E. F. 2006, 

A&A, 453, 163 

j0rgensen, J. K., Johnstone, D., Kirk, H., et al. 2008, ApJ, 683, 822 
Kraus, S., Hofmann, K.-H., Benisty, M., et al. 2008, A&A, 489, 1157 
Lahuis, F., Kessler-Silacci, J. E., Evans, II, N. J., et al. 2006, "c2d 

Spectroscopy Explanatory Supplement" (Pasadena: Spitzer Science 

Center) 

Lahuis, F., van Dishoeck, E. F., Blake, G. A., et al. 2007, ApJ, 665, 
492 

Luhman, K. L. 2007, ApJS, 173, 104 

Manske, V. & Henning, T. 1999, A&A, 349, 907 

Mattioda, A. L., Hudgins, D. M., & Allamandola, L. J. 2005, ApJ, 
629, 1188 

Olivcira, I., Menu, B., Pontoppidan, K. M., et al. 2008, ApJ in press, 
arXiv:0810.0829 

Peeters, E., Hony, S., Van Kerckhovcn, C, ct al. 2002, A&A, 390, 
1089 

Pontoppidan, K. M., Boogert, A. C. A., Fraser, H. J., et al. 2008, ApJ, 
678, 1005 

Pontoppidan, K. M., Dullemond, C. P., van Dishoeck, E. F., et al. 

2005, ApJ, 622, 463 
Pontoppidan, K. M., Fraser, H. J., Dartois, E., et al. 2003, A&A, 408, 

981 

Rcsslcr, M. E. & Barsony, M. 2003, ApJ, 584, 832 
Siebenmorgen, R., Kruegel, E., & Mathis, J. S. 1992, A&A, 266, 501 
Terada, H., Tokunaga, A. T., Kobayashi, N., et al. 2007, ApJ, 667, 
303 

Thi, W.-F., van Dishoeck, E. F., Dartois, E., ct al. 2006, A&A, 449, 
251 

Ulrich, R. K. 1976, ApJ, 210, 377 

van Dishoeck, E. F. & Black, J. H. 1982, ApJ, 258, 533 
van Dishoeck, E. F., Dartois, E., Pontoppidan, K. M., et al. 2003, The 
Messenger, 113, 49 



van Dishoeck, E. F. & van dcr Tak, F. F. S. 2000, in IAU 
Symposium (San Francisco: ASP), Vol. 197, From Molecular 
Clouds to Planetary Systems, ed. Y. C. Minn & E. F. van Dishoeck, 
97 

van Kcmpcn, T. A. 2008, PhD in Astronomy, Leiden Observatory 
- Leiden University, P.O. Box 9513, 2300 RA Leiden, The 
Netherlands 

van Kempen, T. A., van Dishoeck, E. F., Salter, D. M., et al. 2008, 
A&A submitted 

Van Kerckhoven, C, Hony, S., Peeters, E., ct al. 2000, A&A, 357, 
1013 

Verstraete, L., Puget, J. L., Falgarone, E., ct al. 1996, A&A, 315, L337 
Visser, R., Geers, V. C, Dullemond, C. P., ct al. 2007, A&A, 466, 229 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003, 
ApJ, 591, 1049 



